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The glueball spectrum within the Hamiltonian formulation 
of lattice gauge theory (without fermions) is calculated for 
the gauge group SU(2) and for two spatial dimensions. 

The Hilbert space of gauge-invariant functions of the gauge 
field is generated by its parallel-transporters on closed paths 
along the links of the spatial lattice. The coupled cluster 
method is used to determine the spectrum of the Kogut- 
Susskind Hamiltonian in a truncated basis. The quality of 
the description is studied by computing results from vari¬ 
ous truncations, lattice regularisations and with an improved 
Hamiltonian. 

We find consistency for the mass ratio predictions within a 
scaling region where we obtain good agreement with standard 
lattice Monte Carlo results. 


I. INTRODUCTION AND SURVEY 

It is the purpose of this paper to present the results of an 
attempt to compute the spectrum of a lattice gauge field 
theory within the Hamiltonian formulation. Our compu¬ 
tational framework is the coupled cluster method which 
has been presented and discussed in detail in Ref. jl^ . 
We have applied our method to the lattice version of the 
SU{2) Yang-Mills theory in two spatial dimensions which 
is the simplest non-trivial non-abelian lattice gauge field 
theory. We consider this model whose dynamics is given 
by a special Kogut-Susskind Hamiltonian as an impor¬ 
tant test case for controlling the coupled cluster method. 
Although there is no experiment to compare with, there 
exist for this model very reliable results for the ground 
state energy (see Ref. @,^) and for glueball masses (see 
Ref. ) which allow a critical test of the coupled cluster 
predictions. 

In the past, several attempts have been made to com¬ 
pute glueball masses within the same model. We men¬ 
tion the light-front appoach of van de Sande et al. (see 
Ref. 1^) and the maximal tree approach of Bishop et 
al. (see Ref. [pd|). 

Calculations similar to ours, i.e. treating the Kogut- 
Susskind Hamiltonian within the coupled cluster method, 
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were undertaken by Llewellyn Smith and Watson (see 
Ref. [^) and Guo et al. (see Ref. 

Within the present work we are able to go far beyond the 
limitations of this earlier work: 

• When expanding the characteristic correlation 
functions within a localized loop basis (see 
Ref. j^), we can push the expansion up to the 
7th order. 

• We are able to check the validity of the compu¬ 
tation of the characteristic coupled cluster matrix 
elements by comparing two indepedent formula¬ 
tions of the local loop space basis expansion (see 

Ref. [0|@). 

• We can compare the predictions of different lat¬ 
tice regularizations which exist in two spatial dime- 
nions: square, hexagonal and triangular lattice. 

• We compute glueballs with all possible angular mo¬ 

mentum guantum numbers allowed within the lat¬ 
tice formulation, i.e. = 0+,0“,2+,2“, 1^ for 

the square lattice and = 0+, 0“, 3+, 3“, 1^, 2^ 

for the hexagonal or triangular lattice. 

• We also test the quality of our results by comparing 
to those with an improved Hamiltonian. 

Taking dimensionless mass ratios we observe a good scal¬ 
ing window (around (3 = 3 for the non-improved case). 
Judging the stability of the results within the scaling win¬ 
dow by testing the convergence with the order of the cou¬ 
pled cluster expansion, by comparing the predictions of 
the different lattice regularisations and the inclusion of 
improvement, we are able to predict the spectrum of the 
energetically low-lying glueballs within certain errors (see 
Fig. 0 ). Within these errors, our prediction of the glue¬ 
ball spectrum agrees with that of Teper (see Ref. [[l^). 
Since our results are very encouraging it is, of course, 
desirable to extend these computations to more phys¬ 
ical models. The SU{2) lattice Yang-Mills theory for 
ZJ = 3-1-1 dimensions has been worked out (see Ref. p0|), 
but the results are too preliminary to be presented. The 
extension of our model to the SU{3) gauge group is un¬ 
der way (see Ref. |^). The inclusion of quarks, i.e. the 
treatment of full QCD, is possible in principle, but tech¬ 
nically difficult. First attempts using the coupled cluster 
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method for such models exist, but only for the Schwinger 
model (i.e. the 1 + 1 dimensional 17(1) case) (see Ref. p). 
Our paper is organized as follows: 

Section || defines our computational scheme, especially 
details of the coupled cluster expansion. 

Section 01 describes the structure of the two local 
loop space expansions (algebraic versus Clebsch-Gordan 
methods) and how the computation of the coupled clus¬ 
ter matrix elements is done. 

Section fV gives some details of how to project localized 
wave function on good lattice momentum and angular 
momentum (only lattice momentum zero is considered 
within this paper). 

Section ^ describes some calculational details, especially 
the truncation procedure which is necessary to fix the 
coupled cluster method. 

Section Vd presents results for the vacuum state. It is 
shown that within the scaling window the vacuum energy 
agrees with the value obtained by the Greens function 
Monte Garlo method (see Ref. |^) better than within 
0.1% showing that the coupled cluster method is thus 
reliable. _ 

Section VR gives the results for the glueball masses. The 
stability of the results within the scaling window is shown 
by testing stability with respect to the or der of the cou¬ 
pled cluster expansion (Subs ection VII A ), with respect 
to improvement (Subsection VIIB) and with respect to 
the different regularizations (Subsection |VII G| ). 

In the Appendix we discuss the different scales in 
the Hamiltonian and the Lagrangian formulation (see 
Ref. 1^ ) and study the possibility to have a scaling region 
around /? = 3 for our framework. 


II. COMPUTATIONAL SCHEME 

Within the Hamiltonian formulation, the computational 
framework of any lattice Yang-Mills theory is given by the 
Kogut Susskind theory which has been discussed in detail 
in Ref. ||^. In summary, the structure is the following: 
The Kogut Susskind wave functions '^{U) = 

I'({7i,..., Un) depend on link variables Ui {I = 1, ..., N) 
which are elements of the gauge group SU{Nc). N is 
the number of oriented links in a lattice (with d spatial 
dimensions). N is related to the hnite volume chosen. 
Our many-body technique allows us to perform the final 
calculations in the infinite volume limit N —> oo. Our 
concrete calculations are done for the case Nc = 2 and 
d = 2. 

The idea of the coupled cluster method is to refor¬ 
mulate the eigenvalue problem = E'i> as equa¬ 

tions for the ground state correlation function S and for 
the excitation operators F, following from the ansatzes 
I'o(?7) = expS{U) for the ground state and 'I'(17) = 
F{U) exp S{U) for excited states. 

If i7 = l/(2a)(X)/c 9^^fc is the Kogut Susskind 

Hamiltonian (we use the notations (Ic) for link-colour 


quantum numbers, Fic for the momentum operators con¬ 
jugate to Ui and V for the plaquette term) this yields the 
non-linear equation 

y^{Sicic + Si,Si,)-^V = ^Fo ( 1 ) 

^ 9^ 9^ 

for S and the linear equation 

y (Fieic + 2SM = ^(F - Fo)F (2) 

Ic ^ 

for the excitation operators F. We use the notation that 
fic abbreviates “link variable derivatives” of functions 

f{uy. 


flc=[Elc,f] 
flclc = [Elc, [Elc, /]] . 


This (rigorous) form of the eigenvalue problem manifestly 
guarantees the correct volume dependencies of both the 
ground state energy Fq and the excitation energies F—Fq 
(see Ref. [p^). 

Within this paper we will consider excitations corre¬ 
sponding to glueballs with arbitrary lattice angular mo¬ 
mentum and parity, but with momentum zero. Thus the 
eigenvalues are interpreted as glueball masses. We may 
then write 


p _ pk,v 

^ ^(7 / j int 

V 

s = n°5int 

where H^^ is the projection operator on states with 
vanishing (lattice) momentum and angular momentum- 
parity k (e.g. k = = 0+, 0“, 2+, 2“, for the square 

lattice). iy,a are magnetic quantum numbers, for more 
details see Section |r^ We abbreviate Hqq = H*’. 

The essential point of the coupled cluster method is that 
it can be shown that the functions S-mt and Fint are local¬ 
ized (see Ref. |Q) so that it is plausible to expand these 
intrinsic functions with respect to a localized, gauge in¬ 
variant, linked, standard basis in the form 

^i„t(C/) = y^aX“(C/) 

cx. 

E^:iU) = J2Ff'x’^{U) 

a. 

and to define the approximation to the eigenvalue prob¬ 
lem by a truneation of this basis. 

Introducing the constant function via = 1, the “pla¬ 
quette function” by putting HoX^ = 4U and using the 
strong coupling structure xTcic ~ ^aX^ which is true 
for any convenient basis, the coupled cluster equations 
(0) and (H) become equations for the coefficients Sa and 

^ OL 
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,S^+J2c^:o°ik = O+)S^S0 
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0 jy 


= ^iE-Eo)F^-'' 


(3) 


(4) 


Here, the coupled cluster matrix elements C]f^^{k) are 
given by the following prescription: 

Define first numbers c^'l.w which are related to the action 
T{u) (u € Ge) of the lattice Euclidean group Ge (see 
Section IV for details) on the basis x“ and its products 
by 


X^T{u)x*= Y. cil^T{w)x‘^ (5) 

a.,w^GE 

Here only those cases have to be considered where the 
functions x^ and T{u)x^ have a common link variable, 
i.e. only a finite number of possible values for the Eu¬ 
clidean group elements u, w occurs in Eq. 

The coupled cluster matrix elements are then related to 
the matrix elements (b) of the representations of the 
lattice rotations b (u = {b,SL), where a = lattice transla¬ 
tion, see Section by 

C2^J'ik) = l{e^-ep-e^) Y Dty{b-^b')cZL ( 6 ) 

«=(!>,a), 

» = (b',a') 


III. THE LOOP SPACE BASIS 


The basis x“ which we use for our calculations is system¬ 
atically generated from Eq. i) in two steps. 

In a first step, certain subspaces of dimension ms,k 
are generated from geometrically inequivalent plaquette 
products. The second step then is the construction and 
the handling of orthogonal basis systems of these sub¬ 
spaces, especially the computation of the coupled cluster 
matrix elements. This is done in two alternative ways, 
by algebraic and by Clebsch Gordan techniques. 

The subspace is generated from a function Aq^ 

which is given by a linked, standard (5-fold plaquette 
product of the type 

= x^T{u 2 { 6 , k))x^... T{us{S, k))x^ (7) 

k = I,.. .,ns 


where u\{S, k) {X = 2,... ,6) are suitable elements of the 
lattice Eucildean group. Eor 5 = 0,1 we define uq = ni = 

~ A I ^^G ~ A, ■ 


The “ordeE of the functions in the subspace jg (j 0 _ 
fined to be S. In general, the quantum number k charac¬ 
terizes the different, geometrically independent possibili¬ 
ties for the construction of the plaquette products of the 
same order whose total number is ns- 
Up to 7th order we have for our two-dimensional SU(2) 
case 

ns = 1,1, 2,4,12,35,129,495 for the square lattice, 
ns = 1,1, 2, 5,15, 53, 235,1125 for the hexagonal lattice 
and 

ns = 1,1, 2,3, 8,17, 54,162 for the triangular lattice, for 
(5 = 0,...,7. 


The functions Aq^ are characterized by simple loop pat¬ 
terns exemplified up to 4th order in Tabels li and p| . 
Here, also the dimensions ms^k of the corresponding sub¬ 
spaces are given. 

Each function ((5, k fixed) generates an orthogonal 
basis r = 1 ,.., ms,k 7 of and taking all quan¬ 

tum numbers (i5,/c,r) and the limit 6 oo, this basis 
becomes complete for the expansion of localized func¬ 
tions. 

There are two ways to define the basis of explicitely. 


A. The algebraic method 

The idea of this procedure is to act on Aq^ with a suitable 
set of commuting Casimir operators of a lattice gauge 
group. Hereby, this lattice gauge group SU{Nc))^ is 
finite dimensional, where M{5,k) is the number of link 
variables contained in 

This yields a generating system for the space ^nd 
diagonalizing the corresponding Casimir matrices defines 
the orthogonal basis. The relation to the generating sys¬ 
tem allows then to compute also the coefficients (see 
Ref. Q for details and examples). 

The drawback of this method to construct the basis 
X“ = x^'^A’’') is that one has to take care of possible lin¬ 
ear dependencies between the generated loop space func¬ 
tions. In previous investigations (see Ref. |]l^ , 1-10) 
this problem was solved by exploiting the Cayley Hamil¬ 
ton relationship between matrices. We have used in 
this connection a much simpler procedure: If a set 
fi{U ),..., fn{U) contains only m {m < n) independent 
functions, the matrix /i(I7^) {i,k = l,...,n) has for 
suitable fixed variables {U ^,..., U”) exactly the rank m. 
Our experience is that statistically chosen variable sets 
{U ^,..., [/") G SU{2) are suitable in this sense. The lin¬ 
ear relation between the functions fi is then easily con¬ 
structed and the dependent functions can be eliminated. 
This works quite well up to the 6 th order. For higher 
orders, however, this method became too unreliable be¬ 
cause of increasing numerical errors when checking the 
linear dependencies in this way. 


3 






IV. THE LATTICE EUCLIDEAN GROUP 


B. The Clebsch-Gordan method 

Each basis function generated by the algebraic 
method can be characterized by a set of angular momen¬ 
tum quantum numbers adjoined to a certain plaquette 
pattern related to the corresponding function For 

the case of a hexagonal lattice, this pattern consists of 
adjoining a certain angular momentum quantum number 
ji to each link variable occurring in Aq^, the possible 
values of ji being given by simple coupling rules. For 
the square (see Ref. |Q) or the triangular lattice, also 
“intermediate” quantum numbers have to be taken into 
account. This can be systemized by replacing the lattice 
by a certain “net” containing only three-point vertices 
and by defining the plaquette pattern within this net. 
We skip here the description of this net, details can be 
found in Ref. p0| ]. 

The set {ji} of these quantum numbers defines = 
as a well defined function of the link vari¬ 
ables Ui in terms of Clebsch-Gordan coefficients, see 
Ref. For convenience, the set {ji} = 

(ji) j 2 , ■ • ■) will be assumed to be infinite dimensional by 
putting ji = 0 for all links of the net not occurring in the 
plaquette pattern related to 

This makes it possible to work directly with the orthonor¬ 
mal basis and to compute the coupled cluster ma¬ 
trix elements with the help of standard group theoreti¬ 
cal methods. Generalizing the notion {ji} also to non¬ 
standard functions, i.e. putting T{u)x^^‘^ = X^^‘ \ the es¬ 
sential formula allowing to compute the coefficients 
via certain 9j-symbols is given by 

{jf}, {jf (8) 

Gfl 

positively 
oriented I 

3i{v,2) J/%,2) ^l(v,2) / ) (9) 

jliv,3) ) 

where j = -y/ 2 j^PT. 

In Eq. (|^) the (in each case finite) sum runs only over 
those sets {jf} which obey for each I the selection rule 
\ji — jf I < jf < ji + jf and the condition jf yf 0 or 

jf ^ 0 . 

In Eq. (^), we have used the notation 
{l{v,l),l{v,2),l{v,3)} for the three links I adjacent to 
the vertex v. 

For the application of the coupled cluster equation in¬ 
cluding truncation one has to know the order of the func¬ 
tion he. to which subspace jt belongs - the 

plaquette pattern of {j;} does not give this uniquelly. 
This order is, however, easily obtained when construct¬ 
ing the basis x^^‘^ by iterative Glebsch-Gordan coupling, 
see Ref. [Q for examples. 


Any explicit calculation involves the treatmant of the lat¬ 
tice Euclidean group Ge which consist of translations 
and rotatons, Ge = Gt Gs Gr, leaving the (infinite) lat¬ 
tice invariant. 


A. Lattice Translations 

Translations x —> x -|- a are given by translation vectors 


a = OiGi -I- 0202 


where ei = ^ q ^ • For Ibe other unit vector we have that 

TT 

3 

TT 

3 

the hexagonal or triangular lattice, respectively. 

All distances will be measured in units of the lattice spac¬ 
ing. The numbers oi, 02 are then integers which are arbi¬ 
trary for the square or the triangular lattice and obey a 
simple restriction for the hexagonal lattice (see Ref. ). 


02 = 


for a square lattice and 02 = 


cos( 

sin( 



B. Lattice Rotations 


General Euclidean transformation are given by a pair 
( 6 , a) mapping x —> 6 x -|- a where b G Gr is a lattice 
rotation. 

For the square lattice, the “point group” Gr has the eight 
elements (bo, 67 ) which may be written in terms of the 
following 2 x 2 matrices 


bn — 

bn-\-4 — 


cos(n|) sin(n|) 
— sin(n^) cos(n^) 

0 I ' 

I 0 


, n = 0,..., 3 . 


For the hexagonal or triangular lattice the group Gr has 
the 12 elements (bo,bn) which we write in the form 


bn — 

bn+6 = 


cos(n|) sin(n|) 
— sin(n|) cos(n^) 

1 0 
0 -1 


bn , n = 0,..., 5 . 


C. Irreducible Representations 

We now list the irreducible representaions (irreps) of Gr 
(see e.g. Ref. |I|). 

For the square lattice (hexagonal or triangular lat¬ 
tice) we have 5 (6) inequivalent irreducible represen¬ 
tations (irreps) 4 one-dimensional irreps labelled 
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Ai, A 2 , Bi, B 2 and one {two) two-dimensional irreps la¬ 
belled E (El, £’ 2 ). In the continuum limit, these irreps 
will correspond to the lowest spins O’*", 0“, 2+, 2“ and 
(0+, and 2^), respectively. Note that, 

in the continuum limit, the eigenvalues for the quantum 
numbers 7+ and J~ have to be degenerate for J ^ 0 
for a rotational invariant Hamitonian. Explicitly, the ir¬ 
reps read for the one-dimensional cases (n = 0,..., 7) 
for the square lattice {n=0, ..., 11 for the hexagonal or 
triangular lattice) 

D^Hbn) = l 
D^'^ibn) = det(6„) 

D^^ibr,) = (-l)"det(6„) 

D^^ibr,) = (-1)" 

The two-dimensional irreps read for the square case 

D^{bn) =bn (n = 0,... ,7), 

and for the hexagonal or triangular case 

D^^{bn)=bn 

=(-l)"&n (n = 0,...,ll) 


For our concrete calculations we restricted the computa¬ 
tion to eigenstates with momentum zero. This simplifies 
the work because IIq commutes with the operators 11 
A basis for the expansion of the correlation functions S 
or F is then generated from the localized basis x“ via 


x“ |a e IV, 0 < :. < dim . (12) 


The states of Eq. (^^ form an orthogonal basis if the 
functions x°‘ orthogonal and if - for the given k - 
states are left out which are projected to zero under II^q. 
This may occur if the loop space pattern related to x“ 
is too symmetric, i.e. if this state is invariant under a 
sub grou p of the lattice rotation group. Tables 0. 0 


and VI give a survey on the remaining dimension of the 
relevant subspaces up to the 7th order. Note that for 
higher angular momenta only calculations in orders with 
at least (5 > 5 make sense. 


D. Projections on Lattice Momenta and Angular 
Momenta 

The group Ge acts on the functions 

T{u)x^^‘^ = 

I' = u-H 

for u G Ge- 

This allows to define projection operators on good mo¬ 
mentum p within the Brillouin zone by 

a 

States with fixed lattice angular momentum may be con¬ 
structed by the operators 

nt = iVfc^E^^(6-i)r(6„,o), (10) 

n 

where the normalization factor = ddmD^/\G e\ guar¬ 
antees the projector property 


The index a in Eq. & corresponds to a “magnetic” 
quantum number since we have 

T{b)ut,^j2^tAmL'- ( 11 ) 

(T' 

Thus for the eigenvalue problem of a rotationally invari¬ 
ant Hamiltonian a can be fixed, e.g. to cr = 0. 

Eq. ( |ll| ) also shows that the operator H^^ generates the 
representation twice if it is two-dimensional. Thus 
the index e is not a “good” quantum number, the Hamil¬ 
tonian will in general mix states will different values of 

V. 


V. COMPUTATIONAL DETAILS 

With the definition of Eq. (1) for the basis expansion of 
the correlation functions S and E, it is straightforward 
to show the validity of the form Eq. of the coupled 
cluster matrix elements. 

For the numerical study, we have set up a computer pro¬ 
gram which organized efficiently the handling of the loop 
space and of the lattice Euclidean group. For the com¬ 
putation of the coupled cluster matrix element, we used 
up to 6th order both the algebraic and the Clebsch Gor- 
dan method. We found a perfect agreement which was a 
comforting assurance of the validity of our computational 
methods. Details are described in Ref. [p0| . 

For the concrete calculations, we still have to define a 
truncation prescription. For the results presented below, 
we follow the proposal of Guo et al. (see Ref. [^|-Q) which 
yields actually the simplest set of equations. In this case 
one puts in the order S 

cZ:Z=0 for <5(x“^) + 5(x“^)><5. (13) 

Since the expansion of the function S starts with <5 = 1, 
our calculations go up to 8th order if the expansion of 
the correlation functions is pushed to 7th order, as we do 
within this work. 

We have also studied alternative truncations (see 
Ref. p0|), but we could not find any prescription which 
was better than Eq. dil)- 

VI. RESULTS FOR THE GROUND STATE 

It is a special feature of the Hamiltonian formulation that 
it provides both the energy and the wave function of the 
vacuum state which is the ground state of the Hamilto¬ 
nian. 
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Standard lattice Monte Carlo calculations do not give 
results for the vacuum energy density. There exist, how¬ 
ever, computations within the strong coupling expansion 
(see Ref. 0 ) and very reliable Green’s function Monte 
Carlo results (see Ref. !)■ 

In Figure 1^ we compare our coupled cluster results up to 
8 th order to the results of the other methods. As coupling 
variable we use the standard expression j3 = 4/^^. We see 
that we obtain excellent results (better than 0.1%) in our 
scaling region {j3 « 3, see below), we have a good quality 
and convergence up to /3 ~ 4, but our method breaks 
down for large f3 where the Green’s function Monte Carlo 
is still valid. 

An important feature of our results is that the validity of 
the coupled cluster method clearly goes beyond the range 
of the strong coupling expansion which breaks down at 
/3 « 3. (Our figure gives the result of the 18th order 
of strong coupling perturbation theory!) We consider 
this improvement as a necessary condition for obtaining 
continuum limit physics. 

We should stress that the determination of the vacuum 
operator S{U) - by solving Eq. i) iteratively - turns 
out to be much simpler and faster than that of the ex¬ 
citation operator F{U). Thus we do not see any special 
difficulty with the fact that the Hamiltonian formulation 
also involves the determination of the vacuum state. The 
coupled cluster formulation deals with this problem ap¬ 
parently quite effectively. 

In this sense our framework seems to be quite orthogo¬ 
nal to the light front formulation which is based on the 
assumption that the simplicity of the vacuum could be 
helpful for the determination of the spectrum. 


VII. THE GLUEBALL SPECTRUM 

The eigenvalues of the solutions of Eq. (^) have the inter- 
pretion of (approximate) glueball masses with angular- 
momentum = k. 


In Fig. 1^ (upper figure) we present mass ratios of the 
different lowest states of the square lattice relative to 
the mass of the 0+ state in 8th order. As an example of 
the convergence with the order we give in Fig. ^ for the 
three highest orders the mass ratios of the 2“ state versus 
the 0+ state. We compare our data with the scaling 
window results of Teper (see Ref. |^]). As an example 
of mass ratios of higher excitations we show in Fig. || that 
of the O’*'* and the 2+* state, because here a comparison 
with Teper is also possible. 

We assume the scaling window to be at (3 « 3. This 
allows us to estimate the glueball spectrum within certain 
errors by comparing the pred iction of the 7th and 8th 
order at /3 = 3 , see Table VIl. 


Also included are the p redicti ons of the improved Hamil¬ 
tonian (see Sub section VHB) and of alternative lattices 
(see Subsection VHC ). 

For the 0“ sector we have to remark that because the 
glueball in question lies in the continuum of two (low¬ 
est) 0+ glueball states the results are not very reliable. 
This can be seen having a look on the wave functions of 
this states. In the 0“ sector extending the calculations 
from 6th to 7th and 8th order a new state occurs (see 
Fig. ^). In the strong coupling limit (/3 = 0) it has an 
eigenvalue of twice the eigenvalue of the 0+ glueball and 
its wave function is exactly the product of two glueball 
excitation operators: the glueball excitation F is given 
by the product of two plaquettes {two-cluster state): 


7 ^0+ (/3 = 0) = □ , 7th order 

Fq- (/3 = 0) = q n , 7th order . 

As is the strong coupling limit, the scaling region (/3 = 3) 
is dominated by two-cluster states: 


Fo (D 3.0) : ^ 25% + 10 67% 10.57% 5.86% 

+ 4&„ + E% + llTw 


A. Results on the square lattice 

An example of the convergence of the lowest O'*" glue¬ 
ball mass up to the 8th order of the coupled cluster ex¬ 
pansion is given in Fig. ^ (upper figure). We see that 
there is perfect agreement with Hamer’s results of the 
strong coupling expansion (see Ref. 0) up to /3 « 2 
which breaks down for larger /3. Also Hamer’s ESGE 
method (see Ref. 0) becomes apparently unreliable for 
/3 > 3. 

Obviously, we are able to estimate in the region (3 k, i 
the value of the glueball mass within errors. 

The problem is whether this region of couplings contains 
already continuum physics. This we judge from the oc¬ 
curence of an approximate scaling window for mass ra¬ 
tios. 


This feature of the wave functions appears in the whole 
0 “ sector and also in the 1^ sector for the 1^* state. 
The quality of the result on the square lattice is now 
tested by comparing to other computational schemes 
yielding equivalent descriptions in the continuum limit. 

B. Results with improvement 

The definition of an improved actioiWsee Ref. 0 ) or 
an improved Hamitonian (see Ref. Jl^) is not unique. 
Within our work, we have studied the simplest choice, 
namely a tadpole improvement of the plaquette part of 
the Hamiltonian. This is given by the following prescrip¬ 
tions (see Ref. |0 ). 

Replace 
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where mq is related to the vacuum expectation value of 
the plaquette via 


“o = (^ol^tr(n)IV'o) 

which, in turn, may be computed using the Feynman- 
Hellmann formula 


(0|tr(n)|0) = —(0|iF + Atr(n)|0)|^^„. 

Improvement changes the scales of j3. Our results show 
that the scaling region is shifted from /3 « 3 to /3 « 1.7, 
but the prediction of mass ratios should not change. 
Repeating the computation of the ground state with im¬ 
provement we find similar convergence properties as be¬ 
fore. The plaquette vacuum expectation value turns out 
to be Mo ~ 0.84 for /? « 2. 

As an example for the quality of the convergence with im¬ 
provement we show in Fig. y (lower figure) the behaviour 
of the lowest O'*" glueball mass and in Fig. H (lower figure) 
the structure of the scaling window which appear to be 
improved, indeed. 

Taking the predictions of 7th and 8th order at j3 = 1.7 
we ob tain for the glueball ratios the numbers given in 
Table VII which show a remarkable agreement to the 
non-improved results. 


C. Results for alternative lattices 


We have also computed the ground state and the glueball 
spectrum for the hexagonal and for the triangular lattice. 
Using the “geometrical” rescalings 


/3 

Eo 



we find in both cases perfect agreement of the vacuum 
energy to that for the square lattice. 

Using the same rescaling for the glueball masses, we ob¬ 
serve again an approximate scaling region for /3 « 3. 
Since for these new lattices, the scaling quality for the 
2 + glueball is better tha n th at for the O'*' case, the num¬ 
bers presented in Table VII were estimated of the mass 
ratios of the different states relative to the 2+ glueball 
and than rescaled to the O'*" state. 

For both new lattices, we find a new prediction for the 
3"^ and 3“ glueball which turns out to be approximately 
degenerate in the scaling region, as it should be. 

For the triangular lattice, because of level crossings be¬ 
tween the lowest and the first excited state in the I^ and 
in the 3“ sector we find results which differ from other 
lattices. 


Like on the square lattice also on the hexagonal and 
the triangular lattice the 0“ sector is dominated by two- 
cluster states. In addition the two-cluster states influence 
also the first and second excitations in the 1^ and the 3^ 
sector on the hexagonal lattice. 

Nevertheless we perform a “final” estimate of the glueball 
spectrum in units of the 0"*" glueball by calculating mean 
values of the numbers of different orders and different lat¬ 
tices of Table VII, yielding the predictions of Fig. ^ (left 
values) . The errors are deviations of the mean values. 
We do not take into account read off errors. 

We see that, indeed, we obtain predictions quite con¬ 
sistent with the Monte Carlo results of Teper, see 
Ref. [Q (middle and right values in Fig. y). 
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APPENDIX 

Within this Appendix we comment on the question 
whether it is possible to have a scaling region /3 « 3 
if this region starts only with /3 ^ 4 within the standard 
lattice Monte Carlo method (see Ref. Q). 

The point is that a comparison to our Hamiltonian results 
involves a rescaling of the standard “Euclidean” coupling 
gE relative to the Hamiltonian coupling g. For g —> 0 this 
is given by (see Ref. [||) 

/?£ =/3+.077 + O(g2). 

There is also rescaling of the Euclidean masses Me rel¬ 
ative to the Hamiltonian masses M (“velocity of light 
correction”) 

Me = {1 + 0.0845^ + 0(/))M. 

However, one cannot expect that these formulae hold for 
/3 ~ 3 because the difference between [3e and /3 becomes 
very large in the strong coupling regime itself. The strong 
coupling expansions are completely different in the Eu¬ 
clidean and Hamiltonian formalisms: the Euclidean ex¬ 
pansion has a logarithmic singularity for f^E = 0 which 
is not present in the Hamiltonian framework. 

The rigorous relation between Pe and (3 is not known. 
Consequently for j3 ~ i a, direct comparison of masses to 
Teper’s results appears not to be possible. 
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TABLES 


TABLE 1. Construction of the Loop Space Basis: The 
functions A^*^ generating m^^fc-dimensional subspaces are 
characterized by simple loop patterns. Altogether we have 
84 basis elements in the order 5 < 4 for the square lattice. 


A" 






'^6,k 


A^i = 1 1 

AJ4I = □ 1 

A^^ = 0 1 

= □□ 2 

A^^ = B 1 

A^^ = □□ 2 

= Bn 4 

A^^ = am 3 


A4,1 _ 
I\.G — 

11 

1 

A4,2 _ 

no 

2 

A4.3 _ 

Bn 

7 

A4,4 _ 
I\.G — 

gg 

3 

A^5 _ 

□□□ 

5 


gR 

6 

A4,7 _ 
I\.Q — 

gm 

4 


□ 

nm 

9 

A4,9 _ 

□□ 

□□ 

9 

A4,10 _ 

= Ban 

10 

A4.11 _ 

= RH 

7 

A4,12 _ 

= mnn 

6 


84 




































TABLE II. Construction of the Loop Space Basis: The 
functions generating ma^fc-dimensional subspaces are 

characterized by simple loop patterns. Altogether we have 
94 basis elements in the order 5 < A for the hexagonal lattice. 



ms,k 

\S,k 

'^6,k 

A^i = 1 

1 

l^G — 

1 



a4,2 

■l\-G — 'O' 

2 

> 

11 

1 

A4’3 _ O 
^^G — -w' 

6 



~ -O 

5 

l^G — - 

1 

A^’^- 

^^G — w - 

3 

i\g — 

2 

A4.6 - .T4 

5 



A^’'^ - 

4 

i\g — 

1 

^'■G “ O- 

4 

A 3.2 

Aq - o-'- 

2 

A4-9 _ > 

^'■G ~ 

6 

a3,3 _ 4 

3 

A^-io _ 0 - 

^'-G ~ 

10 

A^’^ - -- 

i\g — 

3 

^^G — 

7 

> 

II 

J 

. 

3 

A4,12 _ . ^ 

^'■G — O. 

6 



^G — 

8 



A^’14 - - 

4 



A4-15 _ 

I^G — . >- 

6 






94 


TABLE III. Construction of the Loop Space Basis: The 
functions A^^ generating ma^fc-dimensional subspaces are 
characterized by simple loop patterns. Altogether we have 
58 basis elements in the order 5 < 4 for the triangular lattice. 


j^S,k 

ms,k 

3 A 


'^6,k 

AOh = 1 

1 

A^.i = 

A 

1 



Ah2 = 

A/ 

2 

Ahi = A 

1 

Ah3 = 


7 



Ah4 = 


3 

A^h = ^ 

1 

Ah3 = 


6 

A^’^ = A7 

2 

Ah6 = 


9 



Ahr = 


10 

A^h = A 

1 

Ah3 = 

A 

8 

A^’^ = A7 

2 




A3.3 = 3 

4 





58 


TABLE IV. Lattice Angular Momenta: The number of ba¬ 
sis elements x°‘ for fhe different spins on the square lattice. 


order 

5 <0 

5 <1 

(5 <2 

5 <3 

5 <4 

5 <5 

5 <6 

S <7 

x“ 

1 

2 

5 

15 

84 

557 

4942 

47751 


1 

2 

5 

15 

84 

557 

4942 

47751 

n°'x 

0 

0 

0 

1 

31 

377 

4220 

45132 


0 

0 

2 

8 

65 

497 

4759 

47073 

n^'x 

0 

0 

0 

4 

42 

425 

4386 

45790 

ni^x 

0 

0 

0 

8 

79 

874 

8841 

92340 


TABLE V. Lattice Angular Momenta: The number of ba¬ 
sis elements x°‘ for the different spins on the hexagonal lattice. 


order 

5 <0 

5 <1 

(5 <2 

5 <3 

5 <4 

5 <5 

5 <6 

S <7 

X 

1 

2 

5 

17 

94 

677 

6430 

67036 


1 

2 

5 

17 

94 

677 

6430 

67036 

n°"x 

0 

0 

0 

1 

35 

464 

5613 

63947 

tf+x 

0 

0 

0 

4 

43 

512 

5728 

64573 

O'x 

0 

0 

0 

6 

52 

573 

5989 

65854 

n'^^x 

0 

0 

0 

8 

91 

1078 

11691 

130381 


0 

0 

2 

12 

120 

1128 

12008 

130924 


TABLE VI. 
sis elements x' 

Lattice Angular Momenta: The number of ba- 
^ for the different spins on the triangular lattice. 

order 

S <0 

S <1 

5 <2 

S <3 

5 <4 

S <5 

5 <6 

5 <7 

X 

1 

2 

5 

12 

58 

341 

2958 

27885 

n“A 

1 

2 

5 

12 

58 

341 

2958 

27885 

n°'x 

0 

0 

0 

1 

21 

249 

2601 

26857 


0 

0 

0 

1 

22 

250 

2603 

26859 

O'x 

0 

1 

2 

9 

41 

324 

2807 

27734 

n'^x 

0 

0 

0 

7 

56 

562 

5393 

54550 

tf^x 

0 

0 

2 

9 

71 

577 

5538 

54695 
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TABLE VII. The mass ratios of various glueball states 
relative to the 0^ state for the square, the square improved, 
the hexagonal and the triangular lattice in 7th and 8th order 


in the scaling region. 



square 
/3 = 3 

(5 = 7 5 = 8 

square impr. 

13 = 1.7 

5 = 7 5 = 8 

hexagonal 
/3 = 3 

5 = 7 5 = 8 

triangular 
/3 = 3 

5 = 7 5 = 8 

0 + 

1 

1 

1 

1 

1 

1 

1 

1 

0 +* 

- 

1.47 

1.43 

1.41 

1.51 

1.42 

- 

1.76 

0 " 

1.78“ 

2 .00“ 

1.82“ 

2 .01“ 

1.65“ 

2 .00“ 

2.13“ 

2.26“ 

0 "* 

2 .01“ 

2.06“ 

2.06“ 

2.07“ 

2.08“ 

2.15“ 

2.18“ 

2.32“ 

0 "** 

2.30“ 

2 .20“ 

2.34“ 

2.23“ 

2.25“ 

2 .20“ 

2.53“ 

2.53“ 

1 ± 

2.11 

2.10 

2.14 

2.10 

2.17 

2.15 

2.07^ 

2.04*’ 


2.23 

2 .22“ 

2.27 

2.25 

2.26 

2.17“ 

2.29^ 

2 .12*’ 


2.27 

2.23 

2.30 

2.26 

2.26 

2.19“ 

2.39 

2.35 

2 + 

1.56 

1.61 

1.57 

1.60 





2 +* 

1.77 

1.87 

1.78 

1.83 





2 +** 

1.78 

- 

1.82 

- 





2 " 

1.68 

1.66 

1.69 

1.66 





2 "* 

1.77 

— 

1.81 

— 





2 "** 

2.11 

- 

2.08 

- 





2 ± 





1.60 

1.62 

1.66 

1.65 

2 ±. 





1.65 

1.83 

1.97 

1.99 

2 ±.* 





1.65 

1.96 

2.03 

2.05 

3+ 





2.05 

2.05 

2.20 

2.27 

3+* 





2.26 

2.16“ 

2.55 

2.52 

3+** 





2.31 

2.18“ 

2.63 

2.61 

3" 





2.06 

2.08 

1.69*’ 

1.69*’ 

3"’ 





2.27 

2.52“ 

2.09*’ 

1.97*’ 

3"’* 





2.29 

2.190 

2.23 

2.21 


“predominantly 2-cluster states 
level crossing 
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FIGURES 



P 

FIG. 1. Ground state energy density as a function of the 
inverse coupling /3 = -^ up to 8th order. For compari¬ 
son: Strong coupling expansion SCE Q and Green’s function 
Monte-Garlo GFMG method Q. 




P 

FIG. 2. Excited O’^ state energy as a function of the in¬ 
verse coupling /3 = -^ up to 8th order without and with 
improvement. For comparison: Strong coupling expansion 
SCE and Extrapolation of the strong coupling series ESCE 
with Shafer and integrated differential approximants Q. 
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13 

FIG. 3. Mass ratios of the = 0 , = 2'*' , = 2 

and the state relative to the 0 ^ state as a function 

of the inverse coupling /3 = ^ in 8 th order without and with 
improvement. For comparison: Standard Monte Carlo SMC 
results for f = 1 (only scaling window result). ^ is the 
anisotropic factor ^ a^^t) = lattice spacing in space 

(time) direction. 


FIG. 4. Mass ratios of the = 2 state relative to the 
0*^ state as a function of the inverse coupling /3 = -^ in 6 th, 
7th and 8 th order to show the convergence with the order. 
For comparison: Standard Monte Carlo SMC results for 
^ = 1 (only scaling window result). ^ is the anisotropic factor 
^ as(t) = lattice spacing in space (time) direction. 



f3 

FIG. 5. Mass ratios of the first excitations = O'*"* and 
= 2’'"* glueball states in 8 th order. For comparison: Stan¬ 
dard Monte Carlo SMC results for ^ = 1 (only scaling 
window result). ^ is the anisotropic factor ^ “ 3 ( 4 ) = 

lattice spacing in space (time) direction. 
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p 

FIG. 6. Mass ratios of the = 0 state relative to the 0'^’ 
state as a function of the inverse coupling /3 = ^ in 6th, 7th 
and 8th order. This sector is dominated by 2-cluster states. 
For comparison: Standard Monte Carlo SMC results for 
^ = 1 (only scaling window result). ^ is the anisotropic factor 
^ as(t) = lattice spacing in space (time) direction. 



FIG. 7. Final estimate of the glueball spectrum in units 
of the 0*^ glueball. The bars in t he left column are mean 
values of the numbers in Table VII, The shaded area reflects 
the estimated error. For comparison: Standard Monte Carlo 
SMC results for ^ = 1 in the middle column and for ^ = 5 
in the right column, where f is the anisotropic factor ^ 
with agf^t) = lattice spacing in space (time) direction. 
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